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We have investigated the attractive Hubbard model in the low density limit for the 2D square lattice 
using the ladder approximation for the vertex function in a self-consistent, conserving formulation. 
In the parameter region where the on-site attraction is of the order of the bandwidth, we found 
no evidence of a pseudo gap. Further, we have observed that the suppression of the Fermi surface 
known to destroy superconductivity in one and two dimensions, when these systems are treated 
using a non self-consistent theory (Schmitt-Rink, et at, Phys. Rev. Lett. 63, 445 (1989)), does 
not occur when pair-pair interactions are included. However, we do find a quasiparticle lifetime 
that varies linearly with temperature, similar to many experiments. Thus, although this system 
has a Fermi surface, it shows non Fermi liquid type behaviour over a wide temperature range. We 
stress that our work uses thermal Green's functions along the real time axis, and thus allows for 
a more accurate determination of the dynamical properties of a model than theories that require 
extrapolations from the imaginary frequency axis. 
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I. INTRODUCTION: 

The high temperature superconductors show remark- 
able deviations from Fermi-liquid behaviour in their nor- 
mal state above Tc, the superconducting transition tem- 
perature. Although this seems to be experimentally well 
established, no consistent microscopic theoretical expla- 
nation has been found. 

There are a number of major differences with usual 
metallic (BCS-type) superconductors on which we con- 
centrate for a theoretical understanding. The first one 
is the low dimensionality. For example, the normal state 
conductivity mainly takes place in the two-dimensional 
(d-dimensions shall be denoted by dD throughout this 
paper) copper oxide planes. The second difference is the 
extremely short coherence length of the Cooper pairs in 
the superconducting state. These are known to be of 
the order of 20 A (3-4 lattice constants), and therefore 
much smaller than in usual superconductors (~ 1000 A). 
This fact, together with the extremely low ("bad met- 
als" ) quasiparticle density leads to pairs which are barely 
overlapping (~ 10~^ pairs / coherence volume) with each 
other. Such arguments were first stated in detail by Ran- 
deria ^ who investigated conditions for a crossover be- 
tween superconductivity and Bose condensation of pairs 
of electrons. The small overlap is claimed to be related to 
a separation of a temperature T* , at which pairing takes 
place, from the temperature Tc, at which phase coher- 
ence and therefore superconductivity is established. In 
contrast to this, usual superconductors have ^ 10^ pairs 
/ coherence volume which is believed to be the reason 
why pairing and phase coherence take place at the same 
temperature, the mean- field Tc. Another unusual prop- 



erty of the cuprate superconductors is the linear resistiv- 
ity in the normal state of the optimally doped materials. 
An ideal Fermi liquid should show a resistivity caused by 
electron-electron scattering near the Fermi surface which 
varies quadratically with the temperature. On the other 
hand, a crossover to a linear resistivity due to phonons 
is expected to take place at much higher temperatures, 
above the Debye temperature. 

Finally, one last observation that we wish to focus on 
is the presence of a pseudo gap. This feature was first 
observed with NMR by measuring the spin-lattice 
relaxation rate for the Cu nucleus. The relaxation rate, 
rp^rp , which is temperature independent in normal met- 
als, decreases strongly with decreasing T even above Tc, 
in the high-Tc-cuprates. The pseudogap has also been 
measured by optical experiments where it is found in 
the temperature dependence of the scattering time ob- 
tained from a generalized Drude theory for optical con- 
ductivity data [^IH]. The momentum dependence of the 
pseudo gap has been investigated with recent ARPES 
experiments, and it is found to be present mainly along 
the (fc, 0)-direction, consistent with the proposed d-wave 
symmetry of the superconducting order parameter. It 
has also been observed by ^sR [Q and STM [|j measure- 
ments. At present, there is ongoing discussion addressing 
the notion that the gap can be related to pair formation, 
or other possible precursor phenomena of superconduc- 
tivity, which takes place at temperatures above that at 
which macroscopic phase coherence is established. 

The above-mentioned ideas have led us to consider a 
simple model system in 2D which can describe short co- 
herence length pairs which might exist as preformed pairs 
above Tc. This model, the attractive Hubbard model, 
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for on-site, s-wave pairing only, allows us to focus on 
many of the above-mentioned properties. In particular, 
we have examined the dynamical properties of this model 
to see if the attractive Hubbard model in 2D possesses a 
pseudogap. Further, we have examined the temperature 
dependence of the imaginary part of the single-particle 
self energy to learn how the scattering rate behaves. 

We have focussed on this model in the low density 
limit. When studied in the limit of low band filling the 
attractive Hubbard model represents a system with low 
quasiparticle densities, and therefore the weakly overlap- 
ping pairs proposed to characterize the cuprate super- 
conductors are a natural consequence of this problem. 
An approximation which works well in the dilute limit 
is the ladder approximation; this formalism accounts for 
all possible scattering events between particles that can 
occur in the particle-particle channel (only particle-hole 
scatterings are ignored). 

In a simple and elegant paper Schmitt-Rink et at, 
studied such model systems and concluded, at least in 
a non self-consistent treatment of such problems, that in 
2D a stable, 2-particle bound state persists down to T=0, 
and this leads to a T=0 Bose condensation of composite 
bosons (two fermions pair to form a boson). The physics 
of this phenomenon is that the appearance of preformed 
pairs leads to the elimination of the Fermi surface (there 
are no fermions left), and therefore superconductivity is 
suppressed in favour of T=0 Bose condensation. 

The consequence of this work to the model system un- 
der consideration is as follows: in a non self-consistent 
treatment of the attractive Hubbard model at low densi- 
ties, employing the ladder approximation, the system is 
always unstable towards T=0 Bose condensation of pre- 
formed pairs into an infinite lifetime, two-particle bound 
state 1^. One of the focuses of this paper, and a sec- 
ond motivating force behind our study, is to test if this 
idea survives when the theory is solved in a fully self- 
consistent fashion. That is, when one includes interac- 
tions between pairs, does the physics of Ref. survive? 
A careful and detailed study of the interaction between 
pairs was given by Haussmann ||l^,|ll[, and here we will 
consider these ideas as applied to the observed normal 
state anomalies, including the formation of a pseudogap, 
and to the physics of Ref. |^ . 

The numerical work of Haussmann in 3D suggests that 
the physics of this problem is very different than that of 
its non self-consistent counterpart. In particular, he sug- 
gests that the bound state strongly hybridizes with the 
two-particle scattering continuum, and that this greatly 
reduces the appearance of preformed pairs. Such a ten- 
dency has also been suggested by the work of several 
other authors. Fresard et al. [|2| found in 2D, by ap- 
plying the selfconsistent ladder approximation, that in 
the low doping regime the Fermi-liquid properties are 
fully recovered and that only in the strong coupling 
regime can deviations from Fermi-liquid behaviour be 
expected. (The work of these authors differs from our 
own work in that they use a very different method to 



obtain the dynamical properties of this model in this ap- 
proximation.) Also, Singer et al., employed quan- 
tum Monte Carlo methods and concluded that the two- 
particle bound state, which they refer to as a "band of 
pairs" , is strongly overlapping with the one particle con- 
tinuum, and only at very large attractive interactions 
does it become well separated from the one particle con- 
tinuum. Lastly, recent analytical work of Kagan et al. 
iQ shows that when pair-pair interactions are included 
in a T=0 calculation, in the dilute limit the gap between 
the two-particle bound state and the one-particle contin- 
uum starts to close. 

Our results are consistent with the above-state trends, 
and we do not find that the physics of Ref. survives 
when pair-pair interactions are included. Further, we 
find no evidence of a pseudogap. We have used double- 
time Green's functions to study this system at T>0, and 
thus unlike other studies that examined this system using 
imaginary times (the conventional Matsubara frequency 
formulation) who had to rely on Fade approximants or 
maximum entropy techniques, we are able to examine 
the dynamical properties of this system directly. Thus, 
we believe that we have a somewhat more reliable repre- 
sentation of the dynamics of this system. 

We organize our paper as follows. In we intro- 
duce the model and describe different levels of approxi- 
mation which can be used for the ladder approximation. 
In the next section we introduce a k-averaged method 
along with a generalized spectral representation of all 
temperature-dependent Green's functions which enables 
us to obtain results in a fully self-co nsistent calculation 
along the real time axis. In §[IIB we present numeri- 



cal evidence that aids in justifying these approximations. 
Our main results are presented in 
present our conclusions. 



and in we 
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II. MODEL 



The model Hamiltonian which we consider is the at- 
tractive (negative-U) Hubbard model, given by 

H ^-t (c,>Cj,^ + h-c.) - \U\ »^»,T"^i : 

(1) 

where t is the transfer integral between neighbouring lat- 
tice sites cJ^jCi^CT are electronic creation, annihilation 
operators, respectively, and | i7 | is the strength of the 
on-site attractive interaction between two electrons oc- 
cupying the same lattice site. Throughout this paper we 
restrict our attention to d-dimensional hyper-cubic lat- 
tices. 

For completeness, we begin by reviewing the well stud- 
ied ladder approximation to the Bethe-Salpeter equation 
|^,|o|. The Dyson equation which has to be solved to 
obtain the full one-particle Green's function contains a 
large number of complicated diagrams and cannot be 
solved exactly. However, in the ladder approximation, 
which accounts for all possible scattering events between 
particles that can occur in the particle-particle channel 
(only particle- hole scatterings are ignored), one can find 
a solution for the single-particle self energy. The ladder 
approximation can be motivated in the dilute limit by 
taking kp a, the Fermi momentum multiplied with the 
scattering length (which in 3D is given by a = "'j^^f' , 
which is the effective range of the attractive potential) 
as an additional small parameter. This is so because 
all diagrams which include more than one hole propaga- 
tor (crossing diagrams) are neglected, and therefore this 
approximation is valid in the low density limit. The re- 
peated scattering enters the equation through the vertex 
function T^K^iVln), and since the interaction | [/ | is 
constant, the Bethe-Salpeter equation which determines 
r{K,iil.n) becomes exactly solvable. 

To display this solution we introduce the pair suscep- 
tibility, x(K,ir2„), given by 

m.k 

(2) 

where G(k, ium) is the one-particle thermal Green's func- 
tion. Note that the sign of this function is chosen to be 
different by different authors, and readers should take 
note the consequence of this sign choice in future equa- 
tions. Then, we can express the solution for the vertex 
function as 

f(K, zao = - 1 [/ 1 / (1+ 1 (7 1 x(K, , (3) 

which is shown diagrammatically in Fig. ^ 




FIG. 1. Diagrammatic representation of Eq. (g). The ver- 
tex function T contains the repeated scattering of two parti- 
cles. 

In our work we have chosen to subtract the interaction 
strength from F to obtain an analytic function r(K, z) 
with properties appropriate for examination by a gen- 
eralized Kramers-Kroning analysis (this is equivalent to 
subtracting the Hartree potential from the self energy), 
and thus introduce 

F(K,ia„) = f(K,il]„)-(-|C/|) 

l+\U\ x(K,if]„) ■ ^ ^ 

From now on, we shall only refer to F. 

The criterion which determines a breakdown of the 
normal state due to superconducting pair formation with 
decreasing temperature is known as the Thouless crite- 
rion . We have examined this condition in detail, and 
will report on our results elsewhere. Here, we simply re- 
mark that the Thouless condition is associated with the 
occurrence of a two-particle bound state (with infinite 
lifetime) at the chemical potential, and is signified by 

x(K = O,z = 0) = . (5) 

This equation, and the associated normal state prop- 
erties, can be examined in a variety of increasingly more 
accurate approximations, and we review these approxi- 
mations before proceeding to our results. 

A. Non self-consistent, non-conserving theory 

The self energy in this case (from now on denoted by 
NSCNC) is given by 

S°(k,iw„) = ^^F''(k-fq,iw™ + icj„)G°(q,iw„0 . 

m.q 

(6) 

The index ^ indicates the use of free Green's functions. 
The full Green's function in this approximation is: 

G(k,io.„) = G°(k,zc^„) 

+G"(k,zc^„)I]°(k,ic^„)G'°(k,ic^„) , (7) 
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Most importantly, Schmitt-Rink et al, have used 
this level of approximation and have shown that for any 
2D system with an attractive interaction the system is 
unstable to the effective emptying of the Fermi "circle" 
into the two-particle bound state. Then, a condensation 
at T=0 of these noninteracting composite bosons takes 
place. While some criticism of this simple and elegant 
idea exists [0 , we believe that for any non self-consistent 
theory the ideas of Ref. |^ are solid. What happens when 
one includes self consistency, namely when one includes 
pair interactions, is one of the focuses of this paper. Fur- 
ther, unlike Ref. we believe that it is more appropriate 
to assess the credibility of the ideas of Ref. |^ using the 
same kind of formalism. 



B. Non self-consistent, conserving theory 



S(k, iun) = — ^ r(k + q, iujm + iuJn)G{(i, iujm) (9) 

G(k,ic^„) = (G°(k,zc^„)-i-S(k,zc^„))"' . (10) 

Eq. (|l^) is shown diagrammatically in Fig. ^, and from 
now on we refer to this level of approximation as SCC. 
(Note that we find that nothing new is learned when 
one examines the self-consistent, non conserving level of 
approximation, and in this paper we ignore such equa- 
tions.) The significant achievement of a self-consistent 
calculation is the inclusion of pair-pair interactions, as 
was stated previously by Haussmann pO| ; unlike Hauss- 
mann, who studied a 3D continuous system, here we 
study a 2D lattice system, and are thus able to critique 
the physics of Ref. M . 



We are studying a lattice model which has particle- 
hole symmetry. Unfortunately, the NSCNC approxima- 
tion violates this symmetry. It can be restored if instead 
of Eq. (0) for the Green's function we use 

G(k,ic^„) = (G"(k,zt^„)-i -SO(k,iw„))"' (8) 

The diagram for the full Green's function is shown in Fig. 
H and from now on we refer to this level of approximation 
as NSCC. 
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FIG. 2. Diagram for the single-particle Green's func- 
tion (solid line) in the non self-consistent, conserving ap- 
proximation (NSCC). The thin solid lines represent the 
non-interacting Green's function. 

It was suggested by Serene that the inclusion of 
these new diagrams fundamentally changes the physics 
of Ref. At least for the attractive Hubbard model, 
we do not agree with this claim, and instead find that the 
physics of Ref. Q is still correct and, in fact, is greatly 
simplified when one uses a NSCC. 

C. Self-consistent, conserving theory 

One may treat the ladder approximation fully self con- 
sistently by calculating the self energy using the full in- 
teracting Green's function. That is, one solves for 
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FIG. 3. Diagram for the single-particle Green's function in 
the self-consistent, conserving approximation (SCC), with the 
same notation as Figs. 1 and 2. 

Equations (||,^,^,^ have to be solved in an iterative 
way until self consistency is achieved. Our procedure for 
accomplishing this, along the real time axis, is discussed 
in the next section. 



III. CALCULATION PROCEDURE 
A. Formalism: 

In this section we outline how we solved the self- 
consistency problem in our k-averaged approximation. 
Also, we have employed another approximation in the 
spectral representations for certain functions, and here 
we make this approach clear. Lastly, as was men- 
tioned earlier, our work on this problem involves the real- 
frequency formulation of the thermal Green's functions. 
This formulation follows naturally from the analyticity of 
the retarded (and advanced) Green's function which is, 
e.g., explained in the work of Zubarev |l9| , and we refer 
the reader to this reference for further details. The work 
in this section makes clear the relationship between this 
formalism and the Matsubara frequency formalism. 

In order to solve for the non self-consistent version of 
the ladder approximation, in either conserving or non- 
conserving theories, one uses the lattice Green's function 
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to evaluate the pair susceptibility. The susceptibility is 
used to calculate the vertex, from which one can evaluate 
the self energy. In order to solve for these equations self 
consistently, one uses the self energy to evaluate a new 
approximation for the Green's function, and repeats the 
above process until the resulting Green's function con- 
verges. 

In our k-averaged approximation, we require only the 
momentum averaged pair susceptibility; thus, only the 
momentum averaged Green's function is required. That 
is, denoting the k-averaged quantities by overlined quan- 
tities 



N(m) 



x(*r!„) = l^x(K,zf7„) 

K 



(11) 



(12) 



As a computational approximation, we use a spectral 
representation for the k-averaged functions: they are ap- 
proximated by a number of (5-functions (typically sev- 
eral hundred) along the real axis. These i5-functions 
were placed in a fashion such that they are exponentially 
dense around the chemical potential (with a spacing AE 
^ ksT) and around the lower band edge. The advantage 
of this representation is that all frequency summations 
appearing in our k-averaged approximation can be done 
analytically. 

According to our spectral representation, the averaged 
one-particle Green's function reads: 



1 



■duj 



N(m) 

E 



. iuJn - b f 

3 ^ 



(13) 



To express the Green's function entirely by a series of 
poles along the real axis where to every frequency can be- 
long a superposition of different degenerate energy levels, 
is usually called the Lehmann representation [E0U24|. In 
this way we also understand our approximation above. 

The function A{ljj) is the imaginary part of the (re- 
tarded) one-particle Green's function along the real axis 
p9| , p2[ . The superscript (to) in this and future quantities 
labels the number of the iteration step as convergence to 
self consistency is performed. Further, N{rh) is the num- 
ber of 6 functions which were used, 6™ is the position of 
each 6 peak on the real axis, and a™ is the weight of this 
peak. 

If we define the averaged Green's function in this way, 
we can calculate analytically the frequency summation 
which is needed to obtain the averaged pair susceptibility 



N(m) 



N(rh) 



IUJ„ 



" J ^1 



j.k 3 



tanh 



Pbf 



tanh 



fPbl 



which can be abbreviated as: 



M{rh) 

E 



V 2 



(14) 



(15) 



viz., is of an identical formal structure as the k-averaged 
Green's function in our spectral representation but the 
spectral weights of the poles are now given by c™ at posi- 
tions along the real axis. Here M{m) = 
is the new number of poles that are included in the rep- 
resentation for the pair susceptibility. 

With x^{i^n) we can now calculate r'"(iri„), which 
in our theory is given by 



r(ir2„ 



K 



(1+ I (7 I xiin^)) 



x(i+ \u\x) 



(16) 



We assume that the second term in the expansion, pro- 
portional to the mean-squared fluctuations of the pair 
susceptibility, and higher order terms, can be neglected 
based on our knowledge that they should tend to zero in 
infinite spatial dimensions. In order to obtain a spectral 
representation for T we use partial fractions: 



r (zr!„) = 



jj2 j2^^^™-'^ 



M{f, 



M(m) 

E 



(17) 



The position of the poles /i™ along the real axis is given 
by the zeroes of the polynomial 

M{m) 

l[ix~ d^)~u J2 cfY[{x-dr)^o (18) 

k k l^k 

which have to be determined numerically. The weight 
factors 5™ follow from the equation: 



A/(m) 

k l^k 



(19) 
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With the result for F we obtain an equation for the 
spectral representation of E. 



M{m) N{m) 



a; 



M(m) N(rn) 

EE- 



1 



- hf + hf V, 1 + 
Again we can abbreviate this as 

N{fn) M(m) ~ 



(20) 



E (iw„) 



E 



(21) 



To determine the equation for G in the next level of 
iteration we use partial fractions again: 



N(Q) 
(iuJn) = ^ - 



N(a) 



-E 



pole whose weight is smaller than a certain boundary 
I'iuj) which is again smallest at the chemical potential. 
This has to be done in the way that the loss of spectral 
weight is distributed onto all other poles in order to fulfill 
the sum rules. 

In this way we have solved Eq. (|l^ to Eq. ( p^ ) un- 
til a stable self-consistent solution is obtained. N{0) is 
typically chosen to be around 20, and we end up with 
N{mmax) of the order of 300 and at no step of the cal- 
culation does the number of poles exceeds 3000. 

We can therefore do the whole self-consistency loop by 
calculating all quantities along the real axis. In order to 
show that this is indeed equivalent to the imaginary axis 
Matsubara frequency formalism we show how one can 
calculate the particle number in both formalisms. The 
expectation value of the particle number, < n >, can be 
calculated either by summing over the poles along the 
imaginary axis or by summing over the S functions along 
the real axis: 

<7i>=^ G{iuje) = 1 + (25) 



Af(m+1) 

E 



C{iuJn) - D{iuJn) 
-,n+l 



iojn — b. 

Where we used the abbreviations 

C{iU!n) = (it^n - 6°) Y\_i'^^' 



ri+1 



(22) 



C) 



N(m) M{m) 



D{iu}n) 



E 



C) 



(23) 



The poles 6""*"^ for the spectral representation of 



G {iuin) are given by the zeroes of the polynomial 



N(m) M{m) 



i=-oa 



'=0 



which can be expressed by using the spectral representa- 
tion as an integration along the real axis: 



< n > 



lim 



lim 



1 



1 



:G(z) dz 



1 

2TTi f 1 + el^" TT 
1 /■"" A{x) 



A{x) 



dx dz 



1 



dx 



(26) 



This can be reduced by approximating the spectral func- 
tion as a sum of d functions as done above: 



N{m) 

A{x) « Yl 



af TT S{x 



(27) 



tf) 



E/ ]^(it^n — t™') — Together we can express the particle number with these 

d functions: 



(24) 
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and the weight factors a""*^^ are obtained by inserting the 
results for the poles in Eq. (^). 

By going through one loop of self-consistency in this 
equations the number of poles is increased from 7V(m) to 
iV(0)iV(m) M{m) = N{0)N{m)^{N{m) + l)/2. To avoid 
that the number of 5-functions we have to deal with ex- 
ceeds the number we can handle numerically and to avoid 
divergences which can occur if two poles come too close 
to each other we apply two additional approximations. 
First, we unite two delta peaks to a single one if their 
positions come closer to each other then e(cj), where e is 
smallest at the chemical potential and at the lower band 
edge of the unperturbed system. Second, we neglect a 



< n > « - 

TT 



1 



N{f, 



af TT 5{x - bf) dx 



N(m) 

E 



1 I /3 b" 

1 + J 



(28) 



Note that the dimensionality of the system enters into 
this calculation as the shape of G", whose imaginary 
part is the spectral function of the uncorrelated sys- 
tem. To obtain numerical results we composed the unper- 
turbed Green function G° of a number (typically 20) of 
(5- functions. Since the most important physics is happen- 
ing at the chemical potential and at the lower band edge 
we choose the distance between the (5-peaks to be smallest 
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at these points. To make sure that the distance between 
two i5-peaks is always smaller then kBT, at these points 
we made these distances exponentially small by sampling 
the points with an tanh~^(/3 * a;)-function around the 
chemical potential and around the lower band edge. 

B. Justification of the k-averaged method: 

In the following we discuss limiting cases where the k- 
average approximation made in Eq. ( |l6| ) becomes mani- 
festly justifiable. 

This is the case for large temperatures since 1 — 
(^K/2+q) ~ (^K/2-q)j cis the numerator for the calcula- 
tion of X, is small. Therefore x becomes small. Of course, 
this is simply the uncorrelated limit (note that a k de- 
pendence still survives in the Green's function, since the 
free Green's function has such a dependence). The ap- 
proximation also holds true trivially for small bandwidth 
and becomes exact if the band can be approximated by 
a ^-function [viz., the atomic limit, for which t — > 0). 
In this case no k dispersion is present in the problem. 
Further, the case | J7 | = is trivially correct. 

For I U I not too large, F is determined by the pole 
of the bound state at x = ^/ \ U \, which means that a 
weakly dispersive two-particle bound state is well sep- 
arated from the continuum. The dispersion becomes 
weaker for larger | U \ since the effective transfer of a 
pair is t^ff ~ /U. So, in this case the average over 
k-space is also a good approximation, since all we are 
doing is replacing the bound state below the continuum 
by its average. This is shown schematically in Fig. ^. 




FIG. 4. Schematic of the dispersion of the two-particle 
bound state below the non-interacting continuum of the sus- 
ceptibility x(^) ^) (shaded region) as it arises in a NSC for- 
mulation, if the chemical potential is below the continuum. 
The fc-average approximation is related to ignoring the dis- 
persion of this band of states. 



Probably the most interesting limit is that of large spa- 
tial dimensions. It was argued by Metzner et al. |15[| that 
for a system with large dimensions the k dispersion of the 
self energy vanishes, which in the end is very similar to 
what we do. Such an approximation gives several results 
which are seemingly also valid in 2D and 3D; for a review 
see ||. 

The only case where we can carefully critique the qual- 
ity of our fully self-consistent results found using our k- 
averaged method is for a ID model system. Of course, 
as we argued in relation to Eq. (p^, in this dimension 
our results should be least accurate due to large fluctu- 
ation effects. Therefore, such a comparison should be 
viewed as an upper bound of the potential differences, 
and our method should work much better in any higher 
dimension. 

In the ID case we can directly compare results ob- 
tained by our fc-averaged method (we started with a den- 
sity of states A°{uj) = \/{'k\JW^ - lo"^) and 20 initial 
(5-functions) with fully self-consistent calculations (ob- 
tained by summing over the Matsubara frequencies — we 
used 40 to 60 lattice points and 150 fermionic and 299 
bosonic Matsubara frequencies). Of course, for the lat- 
ter calculation, in order to perform this comparison we 
must restrict ourselves to and calculate quantities that 
are averaged over the Brillouin zone. 

In Fig. H we have plotted the particle number as a 
function of the temperature for a fixed chemical poten- 
tial obtained from the fully self consistent calculations. 
This was done for the two entirely different calculational 
procedures. As one can see in this figure, we had to 
use a logarithmic scale to show the differences between 
these two results. One curve is obtained by applying 
the fc-averaged method while the other curve is obtained 
taking into account the full fc-dispersion and summing 
numerically over the Matsubara frequencies. Both meth- 
ods describe the same physics, which is different from 
the results obtained from NSC calculations. That these 
curves agree so well, and the fact that we are able to use 
our fc-averaged method to such low temperatures, is very 
encouraging. 
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FIG. 5. For a fixed chemical potential {fj, — -0.9) the par- 
ticle number is plotted as a function of the temperature. 
All quantities are measured in units of half the bandwidth, 
W = At for ID, and a value oi \ U \= W \s used. The 
k-average method is compared with a numerical summation 
on a finite lattice (60 sites) and a finite number of Matsub- 
ara frequencies (150). Both methods give essentially the same 
result. Note the logarithmic temperature scale. 

We compared fc-averaged susceptibilities for three Mat- 
subara frequencies, ri„ G {0,^,^}. The results in 
Fig. |6| are obtained with the fc-averaged method — note 
that the analyticity oixi^) allows us to evaluate also the 
results for the fc-averaged method for imaginary frequen- 
cies. 



FIG. 6. For three Matsubara frequencies ({0, 221, and 
for the same ID system the susceptibility x is compared as a 
function of temperature. The two different approximations, 
full k-dependence numerical sum over Matsubara frequencies 
and k-averaged method are compared. In the upper graph 
the imaginary part of x is shown and in the lower graph the 
real part. 

The curves in Fig. ^ are obtained by calculating all 
quantities by summing over (intermediate) Matsubara 
frequencies and by considering the full fc dispersion dur- 
ing the entire self-consistency calculation; then, the fc- 
averaged is done at the very end of the calculation, af- 
ter self-consistency had been established. The tempera- 
ture dependence of these susceptibility quantities shows 
an excellent agreement for these two completely differ- 
ent calculation methods. Thus, the validity of both the 
fc-average method, and of the approximation discussed 
in Eq. (|lq), seems to be fine, and our method success- 
fully reproduces the full self-consistent calculation even 
in ID, where the fluctuations in Eq. ( p^ are expected to 
be largest. 

To critique the predictions of dynamical quantities pro- 
duced by the fc-averaged method, in Fig. ^ we compare, 
for two different temperatures, the imaginary parts of 
the fc-averaged self energy obtained in Fig. |^a by using 
the fc-averaged method, and in Fig. ubby applying Fade 
approximants, as explained in Ref . ||25|] , for Matsubara 
frequencies where the full fc-dependence has been consid- 
ered. 



FIG. 7. The k-averaged imaginary part of /m(E(Lj — /i)) for 
a ID system is plotted. In the upper graph the k-averaged 
method is used whereas in the lower graph a numerical sum 
over Matsubara frequencies (150 points) and 60 k-points with 
a k-average applied at the very end of the calculation. In the 
latter case Fade approximates were used to obtain real axis 
results. The difference at ksT = 0.35 is due to the failure of 
the Fade approximate at sparse Matsubara frequencies. The 
minimum of /m(E(aj — /i)) is clearly resolved in both cases 
for low temperatures. 



of states. The band fillings are chosen to be below half 
filling (n=0.5 is half filling), and therefore the absence 
of the van-Hove singularity in the middle of our "band" 
is unimportant. We choose constant particle numbers 
(ri=0.1 and n=0.3) and calculate the chemical potentials 
as a function of the temperature. 



A. Single-particle density of states — no pseudo gap 



In a similar fashion we compare in Fig. ^the results for 
the one particle density of states. Again, the k averag- 
ing for the Matsubara frequency method was been done 
at the end of the full self-consistent calculation, while 
the fc-averaged method uses fc-averaged quantities only 
throughout its approach to self consistency. In both of 
these figures we find good agreement between the two 
methods, showing that the more conventional Matsub- 
ara frequency method results are reproduced by our fc- 
averaged approach to the pair susceptibility. 




FIG. 8. The ID density of states for the two different ap- 
iroximations is plotted for the same temperatures then Fig. 
I Note the local maximum at the chemical potential. 



Encouraged by these successes of our fc-averaged 
method in ID, below we consider the system of greater 
interest, two dimensions. 

IV. RESULTS 



Of prime interest in our study is the single-particle den- 
sity of states, since the presence of a pseudo gap should 
be apparent in this quantity |24| . Fortunately, this is a fc- 
averaged quantity, and so it follows immediately from the 
imaginary part of the self-consistent (fc-averaged) single- 
particle Green's function. 




-3.0 0.0 3.0 



FIG. 9. The resulting self-consistent density of states for 
the quasi two dimensional system with n=0.1 and | 17 |= 2 for 
four different temperatures is shown. Note that the density 
of states develops a maximum at the chemical potential when 
the temperature is decreased. For comparison we have shown 
the density of states from the NSCC approximation for fcsT 
= 0.2. All energies are given in units of [^] . 

Our results for n = 0.1 as a function of temperature 
(and which implicitly include the temperature-dependent 
chemical potential) are shown in Fig. |^. For compari- 
son, we also show the NSCC result for one temperature. 
To aid in the understanding of these results, we show 
the behaviour of the chemical potential as a function of 
temperature in Fig. |l^ — this latter quantity will be 
discussed in considerable detail in the next subsection of 
the paper. 



In this section we discuss our numerical results ob- 
tained using the method described above. We focus on 
a quasi 2D_system by starting the fc-averaged calcula- 
tion with a G^{lu) that corresponds to a constant density 
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FIG. 10. Flow diagram of the quasi 2D system (71=0. 1) for 
the chemical potential /i as a function of the temperature. The 
different levels of approximation are shown. For comparison, 
the line for the non-interacting system (U=0) is given, and 
the Thouless line which determines the NSC Tc is shown. 
The region of the one particle continuum of a non-interacting 
reference system is shown, beginning at = — 1 on the ^ axis. 
The parameters were chosen to be | [/ |= 2 and n = 0.1 and 
all energies are given in units of half the bandwidth [^] . 

The NSCC result (which is very similar to the NSCNC 
result) at ksT = 0.2 shows that the entire single-particle 
band is above the chemical potential (seen from roughly 
UbT = 0.5 to UbT = 2.5), and that below the chemical 
potential one finds the effects of the two-particle bound 
state. We stress that the "gap" between these two fea- 
tures is entirely different from that found by Janko, et al. 
(which they proposed to be a pseudo gap) in that 
they study a 3D system which, more importantly, has 
the chemical potential inside the one-particle band. As 
is seen in our Fig. |l^, for our NSCC calculation at n=0.1 
the chemical potential is below the continuum arising 
from the single-particle states. Thus, the gap feature 
seen in our figure near w = is not related to a pseudo 
gap. Instead, we believe that this behaviour is related to 
the analytical work for the NSC theory described in Ref. 
0- 

Most importantly, the fully SCC result shows no ev- 
idence of a pseudo gap. We, in fact, see the opposite 
of a pseudo gap, wherein there is an increase in spectral 
weight at the chemical potential. Unfortunately, it is 
somewhat difficult to say what happens to the chemical 
potential for this density — so, we now consider higher 
densities to clarify this situation. 

Figures 0^ show the analogous results for n=0.3. 
Here, the NSCC result is similar to the above NSCC data, 
only now the feature from the two-particle bound state 



is much more clearly resolved. Also, the physics in the 
SCC result is more easily seen. These figures show quite 
clearly that (i) the chemical potential is in the band (for 
the SCC formulation, this is true for temperatures below 
about 0.9), and (ii) no evidence of a pseudo gap is found. 
Instead, as above, there is an enhancement of the spectral 
weight near the chemical potential. This finding is in 
disagreement with other SCC results for a d-wave pairing 
potential p7|. 




FIG. 11. Same plot as Fig. H but for a density of n = 0.3. 




FIG. 12. Same as Fig. ^ but for a density of n = 0.3. 



B. Suppression of superconductivity in 2D 

As mentioned in the introduction, a simple explana- 
tion of the (potential) suppression of superconductivity 
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in 2D was given in Ref. [^. The physics can be related 
to the behaviour of chemical potential, namely, when an 
attractive potential of any strength is present in 2D (or 
ID) a bound state is formed, and this state "attracts" 
the chemical potential to it. Thus, the chemical poten- 
tial is found below the bottom of the band, and no Fermi 
surface is left. No Fermi surface means that there is no 
Cooper pairing. 

We note that the formulation of Ref. [| is NSCNC. Al- 
though Serene has claimed that the physics of Ref. 
1^ is destroyed when one uses a conserving theory, we 
disagree with this. In fact, for lattice systems, we have 
been able to formally show that the physics of Ref. ||] 
survives, and in fact becomes more transparent, in a con- 
serving theory. We will discuss these results in a future 
publication. 

In contrast to this, our data for the chemical poten- 
tial in a sec theory shows that the physics of Ref. ||] 
does not survive the inclusion of pair-pair interactions 
in a fully self-consistent theory. This is easiest to see 
in Fig. na, wherein the low temperature extrapolation 
of the chemical potential is into the (interacting) single- 
particle band. Put another way, there is still a Fermi 
surface at low temperatures in a SCC theory. (That is 
not necessarily to say that this is a Fermi liquid — see the 
next subsection.) We note that the reappearance of the 
Fermi surface was also seen in the SCC work of Fresard 
et al. 

To better understand this, we now consider the be- 
haviour of this bound state via the SCC vertex function. 
Examining the vertex function T{Q) probably best ex- 
plains the different physics between the NSCC and the 
SCC formulations: For the NSCC calculation, the infinite 
lifetime bound state is seen as a (numerically broadened) 
delta function in Im{T{iV)) (inset in Fig. |l^). However, 
Im{T{n)) for the SCC calculation does not show a delta 
peak below the continuum — instead, at the chemical 
potential, with decreasing temperature an enhancement 
near fi can be seen. Note that in Fig. ^ for ft below 
the chemical potential there is a non vanishing imagi- 
nary part of r(ri) that is negative. Therefore, the en- 
hancement of /TO(r(ri)) can also be interpreted as a life- 
time broadened remnant of a two-particle bound state. 
So, we stress again that for the SCC calculation there is 
no infinite lifetime bound state, opposite to the NSCC 
calculation. When treating the system self-consistently, 
pair-pair interactions are included which lead to a finite 
lifetime of the pairs and therefore the condensation of the 
fermions into a bound state is hindered. 




-2 2 4 



FIG. 13. Plot of Im{r{Q. - 2^)) for the quasi-2D system. 
The inset shows the result for kT = 0.2 from the NSCC cal- 
culation with an delta peak indicating the infinite lifetime 
bound state. The SCC result is plotted for three temper- 
atures (0.5, 0.2, 0.07) and shows no infinite lifetime bound 
state. All energies are again in units of [^] with | C/ |= 2 
and n=0.1. 



C. Quasi-particle lifetime at the chemical potential 

The above results showed that in a SCC theory the 
Fermi surface survives at low temperatures. The natural 
question is then: Is this a Fermi liquid? 

To help us answer this question, we consider the (real) 
frequency dependence of the Imaginary part of the self 
energy S. In a SCC theory it exhibits a minimum at 
the chemical potential, and this qualitative behaviour is 
indeed similar to a Fermi liquid. However, in Fig. |lj 
we have plotted this quantity, which is related to the in- 
verse quasi-particle lifetime at the chemical potential, as 
a function of temperature. The calculated value seems 
to indicate a linear variation with temperature (the ex- 
trapolated value is simply a function of our numerical 
broadening, and in this calculation has no physical sig- 
nificance) . 
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FIG. 14. The imaginary part of the self energy at the chem- 
ical potential as a function of temperature for the model pa- 
rameters I C/ 1= 2 and n=0.1. The straight line is a guide to 
the eye. 

To understand this we note that standard arguments, 
based on, e.g., Fermi's golden rule, predict the be- 
haviour of a Fermi liquid. However, and we argue, that 
if the dominant scattering is between lifetime broadened 
two-particle bound states and quasiparticles, then such 
a linear T behaviour is indeed expected. Crucial to this 
argument is the approximate temperature independence 
of the lifetime of the two-particle bound state, something 
that our numerics make clear. Further consequences of 
this simple phenomenology are presently being explored. 



V. CONCLUSION 

We have examined the attractive Hubbard model in 
2D using a k-averaged method. Further, we have stud- 
ied the dynamical properties of this model using real time 
axis thermal Green's functions — this allows for an ac- 
curate determination of quantities like the single-particle 
density of states, and the energy dependencies of the self 
energy and vertex function. We have compared our ap- 
proach to the more familiar Matusbara frequency method 
in ID, and have found little difference between the results 
obtained by the two methods. However, the k-averaged 
method allows us to reach much lower temperatures with- 
out enormous computational efforts. 

We have included pair-pair interactions, as formulated 
by Haussmann, in a SCC theory. Our results for such 
a theory lead us to believe that the attractive Hubbard 
model in 2D, for a correlation energy roughly equal to 
the bandwidth, does not have a pseudo gap. Instead, 
there is actually an enhancement of the spectral weight 
near the Fermi surface. The nature of the ground state 
properties is uncertain, since we find a linear temperature 



dependence of the quasi-particle scattering rate over a 
wide temperature range. The reappearance of the Fermi 
surface in a SCC theory in 2D was noted previously by 
Fresard, et al. |J2|, but the linear T behaviour is new, 
and deserves further investigation. 
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